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Abstract 

We investigate a neutrino mass model in which the neutrino data is accounted for by bihnear R- 
parity violating supersymmetry with anomaly mediated supersymmetry breaking. We focus on the 
CERN Large Hadron Collider (LHC) phenomenology, studying the reach of generic supersymmetry 
search channels with leptons, missing energy and jets. A special feature of this model is the 
existence of long lived neutralinos and charginos which decay inside the detector leading to detached 
vertices. We demonstrate that the largest reach is obtained in the displaced vertices channel and 
that practically all of the reasonable parameter space will be covered with an integrated luminosity 
of 10 fb~^. We also compare the displaced vertex reaches of the LHC and Tevatron. 
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I. INTRODUCTION 



The neutrino sector is one of the most exciting sectors of particle physics today. Our 
knowledge of the neutrino parameters has increased tremendously during the past decade. 
Both the atmospheric and the solar mass squared differences and mixing angles are known 
to a good precision However, only an upper bound exists on the so-called 

CHOOZ angle, sin^ ^^13 < 0.04 [6|], and the absolute mass scale for the neutrino masses, 
^ nil, ~ 0.6 eV [3]. Furthermore, there is an ambiguity in the sign of the atmospheric mass 
square difference leading to the possibility of two different hierarchies, normal or inverse, for 
the neutrino masses. 

As the neutrino experiments enter a precision phase the need to explore different neutrino 
mass models increases. The most popular model for neutrino masses, the seesaw mecha- 
nism jsl, beautifully explains the smallness of the neutrino masses as compared to other 
fermion masses. Nevertheless, it is difficult to test this mechanism due to its very high en- 
ergy scale. The right-handed neutrinos introduced in the seesaw model have masses of order 
10^^ GeV and are too heavy to be produced in colliders^. On the other hand, models where 
the origin of neutrino masses is related to TeV scale physics have been proposed 10|, llll and 
these might be tested at the future or even present colliders 12|, [13 . 



In this paper we will consider a TeV-scale mechanism for generating neutrino masses, 
namely through R-parity violating (RPV) supersymmetry (SUSY) 14| . In this scenario the 
neutrino and neutralinos mix givin g ris e to the neutrino masses [15|]- We will be constraining 



ourself to bilinear RPV (BRpV) 16|, ll7|, thus breaking lepton number but not baryon 



number. This model can be viewed as the effective theory of a spontaneously broken R- 
parity symmetry 
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The BRpV neutrino mass model has only a few free parameters and 
therefore is very predictive. Furthermore, in contrast to trilinear RPV neutrino mass models, 
the constraints from LEP on the R-Parity violating couplings are automatically satisfied as 
the couplings are small. 

We will study the specific case of anomaly mediated supersymmetry breaking 
(AMSB) (19I . 20]. In this scenario the contribution to the soft-supersymmetry breaking 



terms from the super-conformal anomaly, which is always present, is assumed to be domi- 
nant. AMSB naturally solves the flavor problem of the minimal supersymmetric standard 
model (MSSM), since the masses of the two first generations of scalars are automatically 
equal and the flavor off-diagonal terms are given in terms of the quark Yukawa couplings. 
However, anomaly mediation in its pure form is not a viable theory since without any 
other soft symmetry breaking terms, tachyonic sleptons are present in the particle spec- 
trum. Therefore, one normally adds an universal scalar mass term. This minimal anomaly 
mediated scenario (mAMSB) is the one we will pursue in an extended form where we also 
add bilinear R-parity violation. 
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It has been shown that mAMSB with BRpV can account for the present neutrino data 21 



Like any SUSY BRpV neutrino mass model it predicts the normal hierarchy, since the 
neutrino masses becomes strongly hierarchical. The hierarchy is induced because only one 



^ However, one might see there traces in the properties of the left sneutrinos 
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neutrino mass is generated at tree-level and the other masses are generated at loop level. 
In general the atmospheric mass squared difference and mixing angle are related to tree 
level physics, whereas the solar mass squared difference and mixing angle are established 
by radiative corrections. The introduction of RPV will render the lightest supersymmetric 
particle (LSP), which in mAMSB is normally the lightest neutralino, unstable. Clearly, this 
fact will be very important for the collider phenomenology. Indeed, the standard signal 
with much missing energy expected from R-parity conserving (RPC) supersymmetry will 
be depleted. 

The RPV couplings giving rise to neutrino masses are also responsible for the decay 
properties of the neutralino. Therefore, an important smoking gun signal for these models 
is the strong connection between neutralino physics and neutrino mixing parameters. Some 
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ratios of the branching ratios of the neutralino are related to neutrino mixing angles [23 , 
In particular, approximately the same number of muon as taus are expected along with a 
W-boson because their ratio is given by tangent of the nearly maximal atmospheric mixing 
angle. By measuring the decay properties of the neutralino, which is likely to be done by 
the LHC, a severe test of this model is possible. 

Due to the smallness of the RPV couplings the neutralino will have a long lifetime, but 
short enough for it to mainly decay within the detectors at the LHC. A distinct feature of the 
mAMSB scenario is the near degeneracy of the lightest neutralino and the lightest chargino, 
causing even the chargino to dominantly decay through RPV couplings. Consequently, also 
the chargino will have a lifetime in the range interesting for colliders such that it will travel 
a macroscopic distance but decay before leaving the inner detector. A very interesting and 
unique signal will therefore be the observation of displaced vertices from the neutralino 
and the chargino. Displaced vertices may also be produced, although only arising from 
neutralino, in the minimal supergravity (mSUGRA) scenario and has been shown to give an 
excellent reach of the model parameters [2^. As the mAMSB has two different long-lived 
sparticles, there is a richer set of possible final states. 

The prospects for collider discovery of RPV, responsible for the neutrino masses and 
mixings, in mSUGRA have been thoroughly studied |23|, l2J, |25|, l26|, |27| . Also, the discovery 
prospects for mAMSB with R-parity conservation at colliders have been analyzed 28|, l29|, 



301 ]. Nevertheless, the collider signals within R-parity violating mAMSB scenario have been 



analyzed so far only considering trilinear R-parity violation 3l|] . In this paper we will study 



the BRpV-mAMSB scenario, requiring that the neutrino masses and mixings are generated 
by the RPV couplings. We will analyze various generic SUSY search channels for the LHC. 
In addition we will also determine the reach in the displaced vertices channel. 

We organize the paper as follows. In section [III we will review the model and the main low 
energy constraints. In section IIIII we outline our choice of final state channels and describe 
our simulation of the signals and backgrounds. In section IIVI we present our results, as well 
as, our conclusions. 
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FIG. 1: Various sparticle masses for tan/3 = 10, ^ > 0, and m3/2 = ^0 TeV as a function of mo- 



II. ANOMALY MEDIATED SUPERSYMMETRY WITH BRPV 



In this section we will give a short review of the model we analyze. For further details 



we refer to 



22 



The R-parity conserving soft terms will be assigned according to the 
minimal anomaly mediated scenario. In addition we will allow for bilinear R-parity violation 
which contributes in total with six more parameters. The R-parity violating couplings are 
restricted to values which are consistent with the neutrino data, which severely constraints 
the available parameter space. 

Let us start by describing the mAMSB model. This model can be parameterized by three 
parameters and a sign, 

m3/2,mo,tan/3,sign(/i) (1) 

all defined at the high energy scale Mqut- The gravitino mass, ^3/2 , is much larger than 
the other masses, as this is the only one generated at tree level. All the soft breaking terms 
are proportional to the gravitino mass in pure AMSB. However, as mentioned above, in 
order to avoid tachyonic slepton masses an extra universal scalar mass mg is added to all 
sfermions and Higgs masses in mAMSB. The ratio between the vacuum expectation values 
of Hu and is as usual denoted by tan/3. Finally, the sign of the Higgs/Higgsino mixing 
parameter, /x, is free. Our conventions are such that /i enters the superpotential as —fiHdHu- 
For the explicit relations between the input parameters in ([T]) and the soft breaking terms 



in mAMSB see e.g. 32 



A typical spectrum for the anomaly mediated scenario is shown in Figure [H Throughout 
most of the viable parameter space the lightest neutralino is the LSP, with only a small 
area having the stau or the tau-sneutrino as the LSP. Whether the stau is heavier or lighter 
than the tau-sneutrino depends on the value of tan/3. For sufficient high values of rriQ 
the lightest neutralino will be the LSP. The gaugino masses are proportional to their beta 
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functions, resulting in the unique relationship for AMSB: Mi : M2 : M3 ~ 3 : 1 : 7. Here 
Mi are the gaugino masses; Mi is the bino mass, M2 the wino mass and M3 the gluino mass. 
This explains the first of two distinct characteristics of the mAMSB spectrum; the near 
degeneracy of the hghtest neutralino (also the LSP is most cases) and the lightest chargino. 
As the wino mass is much lighter than the others it will be almost equal to the masses of 
the lightest states, with the neutral state being slightly lighter. The wino-like nature of the 
lightest neutralino is important as its interactions are stronger and it will be more easily 
produced at colliders. The second characteristic, the near degeneracy of the sleptons, is a 
less robust feature, as it is a consequence of the assumption of a universal extra contribution 
to the scalar masses and thus a feature of the minimal AMSB. 

The invariance under R-parity, defined by Rp = (— jg normally assumed in 
supersymmetric models. This is mainly motivated by two requirements: to obtain a stable 
proton and to get the LSP as a dark matter candidate. However, the first requirement can 
also be obtained by different symmetries forbidding only baryon number violating terms; see, 
e.g., 33|. Introducing lepton number violating terms has the benefit that neutrino masses 
are generated in an intrinsically supersymmetric way. We will introduce bilinear R-parity 
violation in order to produce the observed neutrino masses and mixings. Thus, we add the 
following term 

Wbrpv = eiUHu , 2 = 1, . . . , 3 (2) 

to the MSSM superpotential. In order to acquire agreement with the neutrino data, the 
R-parity violating couplings must satisfies -C /x. For consistency also the soft breaking 
terms, 

BieiLiHu , 2 = 1, . . . , 3 (3) 

are added to the MSSM. Thus, six parameters related to break down of R-parity invariance 
are introduced. In general this will give rise to sneutrino vacuum expectation values (vevs), 
which in turn leads to mixing between neutralinos and neutrinos. The four heaviest states 
from the diagonalization of the 7x7 neutralino-neutrino mass matrix will be almost pure 
neutralino states and we denote these by ; /c = 1, . . . ,4. Moreover, we arrange them in 
order of magnitude of the masses, thus, x? is the lightest neutralino. The chargino-charged 
lepton mass matrix is treated in an analog way and the lightest chargino is denoted by Xi- 

Clearly the introduction of RPV has important consequences for the collider phenomenol- 
ogy as it will render the LSP unstable. As mentioned above, the lightest supersymmetric 
particle is normally the lightest neutralino in mAMSB as is the case in mSUGRA. As we 
allow for R-parity violation, also the areas with stau or sneutrino LSP are viable, these 
particles will decay. Nevertheless, since the stau or sneutrino LSP parameter space regions 
are very small, we will not discuss these areas further although they are properly included 
in our analysis. In the following discussion we will assume that Xi is indeed the LSP. 

As mentioned above, the R-parity conserving parameters are calculated as in the mAMSB 
scenario. Thus, these are fixed at the scale Mqut and renormalization group equation (RGE) 
running is used in order to extract the low energy parameters. We will use the program 
SPheno 3J] , suitably expanded to the case of RPV, for calculating the mass spectrum and 
decay widths. For the RPC parameters we use 2-loop RGB's and we include all 1-loop 
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threshold corrections. The unification scale is defined as the scale where the U{1) and the 
SU{2) coupling constants meet. The RPV couplings are only dealt with at the low energy 
scale. In the case of the bilinear parameters in the superpotential this can be done con- 
sistently without any additional assumption as also the modulus of fi is calculated at the 
electroweak scale and the bilinear parameters form a closed system within the RGE evolu- 



tion [35|. The corresponding soft SUSY breaking parameters also form a closed system 35 
but one has to assume that there are additional contributions at the high scale similar to 
the case of the scalar mass parameters squared to get a consistent picture. 

The parameters of eqs. ((21) and (I3l) are determined with the help of neutrino physics. First 
we trade the Bi by the sneutrino vevs Vi using the tadpole equations. The neutrino masses 
and mixings are best parameterized using the quantities Aj = fiVi + eiVd and = Vijej where 
Vij is the mixing matrix of the tree level neutrino mass matrix 17|] . In case the tree-level 
contribution dominates in the effective neutrino mass matrix, the modulus of A is fixed by 
requiring the correct atmospheric neutrino mass difference squared, the atmospheric neutrino 
mixing angle, and the Chooz angle within the allowed experimental range. In addition, the 
ratios |e|^/|A| and ei/e2 are fixed by requiring the correct solar mass difference and solar 
mixing angle within the experimental range, respectively. 

The very particular near degeneracy of Xi ^^id xt preserved in our model, as the 
R-parity violating couplings have little impact on the sparticle masses. It is important to 
calculate the chargino and neutralino masses very precisely, as the mass splitting between 
these particles, Am^ = m^o — rn^+, is very small and the exact value can have important 
consequences for the chargino lifetime. In our numerical evaluation, the neutrino-neutralino 
and chargino-charged-lepton mass matrices are evaluated to 1-loop order. 



A. Constraints on the model 

Here we will discuss the existing constraints from low energy observables. Besides the 
requirement of agreement with neutrino data, there are also important constraints from the 
LEP data, the rare process b —>■ s'-f and the anomalous magnetic moment of the muon. 

We use the following neutrino constraints, which are the present ones at 90% confidence 
level: 

+ 7.3 X lO^^eV^ < Am4i < +9.0 x lO^^eV^ 

0.25 < sin^ e^ai < 0.37 (4) 
1.5 X IQ-^V^ < \Aml^J < 3.4 x lO'W 

0.36 < sin^ e^tm < 0.64 (5) 



It is not always possible to succeed in generating the observed neutrino masses and mixings, 
and such points will be excluded from our analysis. 

The LEP collider at CERN has already put lower bounds on some of the super symmetric 



sparticle masses. We have implemented the following constraints from LEP |36|: 

TUi > 95 GeV , > 85 GeV , > 79 GeV , (6) 
> 95 GeV , m^o > 42 GeV , niho > 95 GeV . 
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We have checked that also the Tevatron bounds on squark and gluino masses are satisfied 
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The hmit on the Higgs mass is somewhat optimistic as the bound in the MSSM is 
the same as for the SM Higgs throughout most the available parameter space ^ . For this 
reason we exhibits the contour for m^o = 114 GeV in our plots of the collider reach. The 
lower bound on the chargino mass translates almost directly to a lower bound on 7713/2 • For 
tan/? = 10 we must require m3/2 > 30 TeV in order to satisfy the chargino mass bound. 
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FIG. 2: The supersymmetric contribution to the anomalous magnetic moment of the muon for 
tan/3 = 10 and fi > (left), ^ < (right panel). The dark (red) area is theoretically excluded, 
due to either tachyonic particles or the fact that the electroweak symmetry is not broken. The 
light (cyan) area in the left panel is the allowed 3<t range. We also show the contour for (5a^ = 
—4.8 X 10^^^ in the right panel, which corresponds to the Aa lower bound. 

The value of the supersymmetric contribution {6a^) to the magnetic moment of the muon 
is shown in Fig. [2] for tan/3 = 10. The most recent data shows a 3.4cr disagreement with 
the Standard Model (SM) value, having Sa^ = (27.6 ± 8.1) x 10"^° H. As is weU known 
the sign of tracks the sign of the /x-parameter, in theories where the wino and the bino 
mass are positive. Therefore, clearly the case /i < is disfavored. 

The value of BR(fe — > 57) is shown in Figure [3] for tan/? = 10. The current experimental 
value reads, BR(6 S7) = (3.55±0.26) x 10~^, which is in fine agreement with the Standard 
Model expectations. The supersymmetric contribution to this rare process is enhanced by 



tan /3 40| and in general the rate only deviates from the standard model value for large 
tan p. Therefore, the whole region in Fig. [3] lies within the 3a range. Please note that as the 
gluino mass is negative in anomaly mediated supersymmetry one does not find a preference 
for /i > as in the mSUGRA scenario. 



Only for small pseudo scalar mass can a lighter Higgs boson be permitted 38 1 
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For large values of tan/3 (about 40) and /i > the neutrino data cannot be fitted for 
low values of mo. Moreover, for large tan/5 the bound on the lightest neutral Higgs, as 
well as the rate of 6 ^ 37, rules out a large corner of the low mo and low m3/2 parameter 
space. Furthermore, for fi < the constraints from g — 2 oi the muon excludes all of the 
parameter-space at 3cr for tan/5 = 10. As the supersymmetric contribution to the magnetic 
moment of the muon is proportional to tan/5, the disagreement with the measured values 
only becomes worse at large tan (3. Henceforth, we will focus our analysis at a relative low 
value of tan/5, which we fix to be 10, and assume that /x is positive. 



B. Neutrino Parameters 

In this section we study a case solution for neutrino physics within the context of BRpV- 
mAMSB. We concentrate in the following point of mAMSB parameter space, 

m3/2 = 40 TeV, mo = 500 GeV, tan /5 = 10, /i > (7) 

which leads to reasonable values for B{h ^ 57) and 5a^, as can be seen from Figs. [2] and 
131 In the spectrum of this model, the Higgs sector is characterized by a light Higgs mass of 
ruh = 111.4 GeV and a relatively heavy charged Higgs with 'mii± = 834 GeV. The LSP is 
the first neutralino with m^o = 120.85 GeV, followed by a nearly degenerate chargino with 
m^+ = 120.88 GeV. The lightest slepton is the stau with m^ = 458 GeV, and the lightest 
squark is the stop with rrit-^ = 672 GeV. 

The main effect in collider physics of the presence of BRpV is the instability of the 
neutralino. But in addition, in the neutrino sector we gain a mechanism for generating 
masses for the neutrinos, which in turn explain their oscillations. An example solution for 
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FIG. 4: in the £2-63 plane for the left frame and in the A1-A3 plane for the right frame. 

neutrino physics, which we call benchmark 1, is given by the following BRpV parameters, 
ei = -0.0117 , 62 = -0.43 , 63 = -0.246 GeV 

Ai = -0.0467 , A2 = 0.00305 , A3 = 0.0689 GeV^ . (8) 

It predicts the following neutrino observables, 

Aml^^ = 2.4 X 10-3eV^ Am^^i = 8.0 x IQ-^eV'^, 
tan2 e^tm = 1-27 , tan^ ^^oi = 0.49 , tan^ ^^eac = 0.027 . (9) 

These results are calculated from the full one-loop renormalized 7x7 neutralino-neutrino 
mass matrix. 

In order to gain some insight into the problem, we study next some approximations. It 
is known that for small BRpV parameters, the 3x3 effective neutrino mass matrix takes 
the form, 

M^. = AAiA, + B{A,ej + A,e,) + Ce,e,- (10) 

where A receives contributions from tree-level as well as one-loop, and B and C are one- loop 
generated. These three parameters depend only on MSSM masses and couplings and not on 
BRpV parameters. All the dependence on BRpV is in the and Aj. From the 7x7 mass 
matrix in benchmark 1, the corresponding numerical values for the A, B, and C parameters 
of the 3x3 effective mass matrix are, 

A ^ -2.10eV/GeV^ B ^ 0.157 eV/GeV^ C ^ -0.162 eV/GeV^ (11) 

The error we make when we use the approximation in eq. (fTOl) can be estimated with a 
evaluated at the input values given in eq. ([H]), where we defined 

2^ / AmL,-2.35 Y / Am^,i-8.15 Y / sin=^ ^atm - 0.51 V / sin^ gsoi - 0.305 ^ 

^ V 0-95 / V 0.95 J V 0.17 J V 0.075 J 

(12) 
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FIG. 5: Partial iii the e2-(-z plane: atmospheric mass for the left frame, and atmospheric angle 
for the right frame. 



where the central values and Sa deviations were taken from ref . [411 • The atmospheric mass 
difference is given in 10~^eV^, and the solar mass difference in 10~^eV^. Both are defined 
to be positive. The neutrino observables calculated with the effective 3x3 mass matrix give 
the value = 5.7, most of it coming from the solar angle, indicating the kind of error we 
make when we use it. 

For illustrative purposes, we find now the least modified values of BRpV parameters that 
give a good solution for neutrino observables calculated using the diagonalization of the 
effective 3x3 neutrino mass matrix in eq. (fTOj) . We call it benchmark 1': 

ei = -0.0117 , 62 = -0.50 , es = -0.16 GeV , 

= -0.064, A2 = 0.00305, A3 = 0.033 GeV^ . (13) 

The difference between this benchmark 1' and the one in eq. ([8]) indicates us how erred 
would be the determination of BRpV parameters if we do not use the full 7x7 mass matrix 
in the calculation of the neutrino observables. 

In Fig. Hlwe plot which measures the deviation of a given model prediction from the 
experimental measurements. We calculate using the 3x3 effective mass matrix, and we 
use benchmark 1'. In the left frame we vary €2 and €3, keeping all the other parameters con- 
stant as indicated by benchmark 1'. In the right frame we vary Ai and A3. By construction, 
the minimum appear at the values defined by benchmark 1'. The fact that small deviations 
on the BRpV parameters produce very large values of indicate us how sensible are the 
neutrino observables to them. 

To explore in more detail the dependence of each neutrino observable on BRpV param- 
eters we define the quantities xh ^ = 1,4 as the calculated using only the z-th term in 
eq. (fT2!) . As in Fig. |U the xl ^^^e calculated using the 3x3 effective neutrino mass matrix 
approximation. We start with atmospheric parameters in Fig. [51 In the left frame we have 
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2.35 



0.95 



(14) 



which is no more than the first term in eq. ( U21) . associated with the atmospheric mass 
squared difference, plotted as a function of e2 and e^. In the right frame we have the solar 
angle represented by xi, as a function of the same parameters. The behaviour of can be 
easily understood with some approximations which we develop next. 

In normal circumstances the A-term dominates over the other two in eq. (fTO!) because it 
receives contributions at tree-level. Nevertheless, depending on the relative values of the 
and Aj parameters, it is possible for the C-term to dominate the neutrino mass matrix. This 
is the case with the example we are studying. We have obtained approximated solutions for 
eigenvalues and eigenvectors of the 3x3 effective neutrino mass matrix when the C-term is 
much larger than the A and B terms. In this case, the neutrino masses are, 



m2 



'e-x (A X e)p 



A)^ 



A 



(15) 



up to terms of second order, while the lightest neutrino has exactly mi = when the mass 
matrix has the form in eq. ( fTOl) . The eigenvectors are given by the following expressions. 



B + A 



V2 



ex (A X e) 
|e X (A X e)\ 
Axe 



B + A- 



ex (Axe) 
e- A 



|e X (A X e)| 



(16) 



|Axe1 

also up to terms of second order. Note that the eigenvectors are orthogonal and normalized 
up to the order we are working. The matrix Upmns is formed with the eigenvectors in its 
columns. 

Using these approximated expressions, and neglecting A2 and ei, we find. 



tan^ ^atm 



mn 



2\2 



tan^ 6. 



sol 



VZ,2 
V3,3 

V2,l 
V3,l 



A' 



£2 

Ai 



A? 
2 



A^ 



1 + (63/62)^ 



(17) 



With these approximations for the neutrino observables, we can easily understand the dif- 
ferent figures. The atmospheric mass squared difference Aml^^ in eq. (fT7j) indicates that 
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FIG. 6: Partial in the A1-A3 plane: solar mass for the left frame, and solar angle for the right 
frame. 

constant values of this observable are obtained at circumferences in the €2-63 plane, which 
is exactly what we see in the left frame of Fig. [5l At the same time, constant values of the 
atmospheric angle are obtained at straight lines, which is confirmed in the right frame of 
Fig. [H From eq. ( fTTl) we also see that the dependence on Aj of the atmospheric parameters 
is weak, and we do not plot it explicitly. 

In Fig. [6] we concentrate on the solar neutrino parameters, and we study them as a 
function of Ai and A3. From eq. f|T7|l we see that the solar mass squared difference has a 
constant value at ellipses in the A1-A3 plane, and they can be seen in the left frame of Fig. [6l 
with the eccentricity of the ellipse depending on the ratio €3/62. Similarly, in eq. (fTTl) we 
learn that constant values of the solar angle are obtained at straight lines in the A1-A3 plane, 
and they can be seen in the right frame of Fig. [61 The dependence of the solar observables 
on is not weak, but we do not show it here. 

C. R— parity violating decays of SUSY particles 

The R-parity violating interactions are rather feeble since they are related to neutrino 
physics. Therefore, the R-parity violating effects are expected to be small except in processes 
suppressed in the R conserving scenario. In fact, the pair production of SUSY particles is 
nearly the same as in the case of conserved R parity and single production of SUSY particles 
is strongly suppressed. The main manifestation of R-parity violation is the fact that the 
lightest supersymmetric particles decays. 

We evaluated all possible R-parity conserving as well as R-parity violating decays for 
all particles. Notwithstanding, with the exception of Xi Xi R-parity violating 
channels are strongly suppressed as they are proportional to the ratios l^p/l/ip ~ 10~^ 
and |Ap/| det(m^o)| ~ 10~^ and, thus, do not play any role in our analysis. Due to the 
near degeneracy of the lightest neutralino and chargino the decay of the latter through R- 
parity violating couplings is significant (often dominant) and we evaluate all RPV and RFC 
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FIG. 7: The branching ratios of the hghtest chargino (left) and the hghtest neutrahno (right) as a 
function of 7113/2 for tjiq = 800 GeV, tan/3 = 10, and /i > 0. We use the symbol 1 for an electron 
or a muon. By q we mean any quark not being the bottom quark, and by the ud combinations we 
mean any set of u-type quarks and d-type quarks. In general we sum over all states in the decay 
channel. 



decay channels. The decay of Xi a-^d Xi ^-^e calculated using the two-body decay whenever 
possible. Thus, for masses above the ly-mass we compute the decays Xi ~^ ^fW^, whereas 
below My/ we compute the full three-body decay into all possible final states. In the three- 
body decays, we include all intermediates states, such as neutral scalars and pseudo-scalars 
which can be important in some regions of parameter space 2J]. However, the effects from 



the scalar intermediate states, except the Standard Model like Higgs, are less import in our 
AMSB model, as the scalars are fairly heavy. 

The decay of the lightest neutralino in our model takes place through Z, neutral and 
charged scalar and squark states. The decays of the neutralino can be classified as leptonic 
{yll\ semileptonic (^gg', i^gg) as well as invisible [yvv). The possible BRpV chargino decay 
channels, induced through the same intermediate states as for the neutralino, are i^g'g, iqq., 
M and Ivv. Moreover, R-parity conserving chargino decays exist with the most important 
RFC channels being xX ~^ Xi^^ ^ind xX ~^ Xi^~^^- The mass difference Am^, is mainly 
dependent on the value of tan /3 and 7713/2, and is so small that it suppresses the RFC channel 
substantially. For this reason the RFC decays will have the largest branching ratio for small 
m3/2 and large tan/3, since in this case Am^ is the largest. 

To understand better the decay of these particles we plot their branching ratios in Fig- 
ure [71 The chargino decay channels with a branching ratio above 10~^ are plotted and for the 
neutralino we plot all channels above 10~^. We have differentiated between bottom quarks 
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(denoted b) and other quarks (denoted q). Also, we distinguish between a muon or an elec 



tron (commonly denoted by /) and a tau. We note, that as has already been observed |23l. |24 
there is an important connection between the neutralino decay and the neutrino parameters. 
In particular, due to the large atmospheric mixing angle, an almost equal number of muons 
and taus is expected to be produced along with a W. This can clearly be seen from Fig. [71 
as the W-mediated channels exhibit this property. We can also learn from this figure that 
for heavier neutralinos, i.e. larger m3/2, the most important decay channels are fiW, tW, 
and uZ with the uh mode having a sizeable contribution. The three-body decays grow with 
the increase of 1713/2, becoming important at high 7713/2 values. 

The chargino branching ratios also exhibit a number of near equalities. For instance, 
from the right panel of Fig. [7] we can see that BR(x^ — s> iZ) ~ BR(x||'' — > tZ) as well as 
BR(x^ — > ih) ~ BR(x^ rh). Indeed, these equalities are found for all values of 1713/2, 
being a consequence of the fact that the same R- par ity breaking parameters responsible 
for the SUSY decays govern also neutrino physics [42]. Furthermore, the chargino decays 
predominantly into i>W for all values of 7723/2 and the other important decays are iZ, rZ , 
Ih and rh. The mass difference Am^ for tan/5 = 10 is at most 300 MeV and as can be seen 
from Fig. [7], the chargino RFC branching ratios are very small and can be neglected in the 
analysis. 

III. COLLIDER SIGNALS 

We will analyze the LHC discovery potentials for BRpV-mAMSB in various channels, that 
is, we study a myriad of channels, ranging from jets+missing energy to multilepton channels 
in addition to the displaced vertices signal. Throughout this paper we use SPheno |34| to 
generate the particle spectrum and decays which are tabulated in the SLHA format j43| . 



The signal and background generation was carried out with PYTHIA 4J] version 6.409 



adopting the CTEQ5L parton distribution function [45 . 



A. Canonical SUSY final state topologies at LHC 

We considered several canonical supersymmetry signals for the LHC, following what has 
been presented in Refs. 28[ and 46 



1. Inclusive jets and missing transverse momentum (IN) : in this class of events we include 
all events that present jets and missing^. In R-parity conserving scenarios this 



channel is one of the main search modes 28 



2. Zero lepton, jets and missing transverse momentum (0£): the events in this class 
present jets and missing without isolated leptons (e^ , /x^); 

3. One lepton, jets and missing transverse momentum (li): here we consider only events 
presenting jets and missing -^t accompanied by just one isolated lepton; 
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4. Opposite sign lepton pair, jets and missing transverse momentum (OS): the events in 
this class contain jets, missing ^t, and two isolated leptons of opposite charges; 



5. Same sign lepton pair, jets and missing transverse momentum (SS): here we consider 
only events presenting jets and missing accompanied by two isolated leptons of the 
same charge; 

6. Multileptons, jets and missing transverse momentum (M£): we classify in this class 
the events exhibiting jets, missing '^t accompanied by three or more isolated charged 
leptons. 

In our analysis we defined jets through the subroutine PYCELL of PYTHIA with a cone 
size of Ai? = 0.7. Charged leptons (e"*^ or /z^) were considered isolated if the energy deposit 
in a cone of Ai? < 0.3 is smaller than 5 GeV. Furthermore, we smeared the energies, but 
not directions, of all final state particles with a Gaussian error given by b^EjE = 0.7/a/E 
{E in GeV) for hadrons and AE/E = O.IS/V^ for charged leptons. 



We perform our event selection along the same lines of 28|]. Initially, we applied the 
following acceptance cuts: 

CI We required at least two jets in the event with 

j4 > 50 GeV and \rij\ < 3 . (18) 

C2 The transverse sphericity of the event must exceed 0.2, i.e. 

St > 0.2 . (19) 

This requirement reduces efficiently the large background due to the production of two jets 
in the SM. 

C3 A potential source of missing transverse momentum in the background is the mismea- 
surement of jets. Therefore, we imposed the following cut in the azimuthal angle between 
the jets and the missing momentum (A(/?) to reduce the background 

^ < < ^ • (20) 



After applying the acceptance cuts C1-C3 we impose further cuts to each final state 
topology. In order to achieve some optimization of the cuts in different regions of the 
parameter space we used floating cuts E^ which can take the values 200, 300, 400 and 
500 GeV. Considering that our main goal is to evaluate the impact or R-parity violation 
interactions on the searches for SUSY we have not tried to further improve our cuts. Given 
one value of E^, we further required for all topologies that 

> Elj. and i)T > E^ . (21) 

where pl^'^ stand for the transverse momenta of the two hardest jets. These are the only 
additional cuts applied on the IN topology. 
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/Background 


IN 


11 


OS 


SS 


mi 


200 GeV 


261. 


178. 


5.4 


0.40 


1.3 


300 GeV 


24. 


17. 


0.32 


0.013 


0.001 


400 GeV 


4.0 


3.1 


0.015 


0. 


0. 


500 GeV 


0.88 


0.63 


0.003 


0. 


0. 



TABLE I: Total background cross section in fb as a function of for the channels considered 
here. 

For the events in the 0£ class we veto the presence of isolated leptons with 

p^T > 10 GeV and \r]i\ < 2.5 . (22) 

In addition to comply with (I2T]) . li events must present only one isolated lepton satisfying 

p^T > 20 GeV and \r]i\ < 2.5 . (23) 

Since the SM production of W's is a potentially large background we also imposed that the 
transverse mass cut 

"^r(|^r,PT) > 100 GeV . (24) 

For the OS (SS) signal we required the presence of two isolated charged leptons with 
opposite (same) charge after imposing (12T]) . The hardest isolated lepton must have 

p^T > 20 GeV and \r]i\ < 2.5 , (25) 

while the second lepton must satisfy ( 1221) . Multilepton events [Mi) must pass the cuts (l2Tll 
and exhibit three or more isolated leptons with the hardest one satisfying (J25ll and the other 
leptons complying (l22l) . 

The most important SM backgrounds to the canonical SUSY searches are 

• the process with highest cross section at small is the QCD production pp jjX 
where j denotes a jet; 

• tt production that contributes to many final state topologies due to its decay into 
WWbb; 

• associated production of weak bosons and jets which we denote by Wj and Zj; 

• double weak boson production VV with V = Z oi W; 

• production of a single top quark. We did not consider the gluon-VT contribution to 
this reaction because it is not included in PYTHIA. 

We present in Table I the total background cross section after cuts for the final state 
topologies that we analyzed. The main source of background is tt production for all the 
process with Wj and Zj also contributing to the IN background. Moreover, Wj also plays 
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an important role in the li topology. For additional information on the backgrounds see 



Ref. 46 



Comparing our results for the SM backgrounds with the ones in Ref. |28|] we can see 
that they agree within a factor of 2. This is indeed expected since PYTHIA and ISAJET 
have different choices for the event generation, specially in the hadronization procedure. 
This does not pose a serious problem for the LHC experiment since the backgrounds can be 
obtained from actual data. 



B. Displaced Vertices at LHC 

In our BRpV-mAMSB model both Xi and Xi can travel macroscopic distances before 
decaying. Consequently, the long lifetimes of these particles can give rise to a further striking 
signal, that is, the existence of displaced vertices in the events. However, the decay must be 
confined within the inner parts of the detector in order for it to be fully reconstructible. 

The decay lengths in the rest frame of Xi ^"^^ Xi ^ire depicted in Fig. [8] as a function 
of the neutralino (chargino) mass for different choices of parameters. We can see from this 
figure that these SUSY particles typically decay within the inner detector, even when taking 
into account the time dilation 7-factor. Due to the fact that the sparticles production is 
completely dominated by R-parity conserving interactions, there will almost always be two 
long lived sparticles in the reaction chain that may lead to a displaced vertex, either Xi 
or Xi ■ Thus, our analysis will require the presence of two displaced vertices in the event. 
Notwithstanding, there is a small corner of the parameter space where Xi Xi decay 
very fast as we can see in the bottom curves of Fig. [HI This region is characterized by the 
LSP being the stau, and consequently, Xi ^^id xt ^^st decay is due to R-parity conserving 
interactions. Moreover, in this region of the parameter space the stau decay is too fast, at 
least by a factor of 10^ with respect to chargino and neutralinos, washing out the displaced 
vertex signal. 

In order to be able to see the vertex, where the chargino or neutralino decays to the 
secondary particles, we must require that the decay products are such that this vertex can 
be reconstructed. The following decay modes allow us to reconstruct the neutralino decay 

h'i~^i~ with i = e, denoted by ii; 
vqq denoted jj; 
rg'g, called rjj; 
vhh, that we denote by hh\ 
UT^T", called rr; 
rui, called ri. 



vertex 

• X? 

• X? 

• X? 

• X? 

• X? 

• xS 
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(GeV) 

FIG. 8: The decay length in the rest frame of the lightest neutralino and the lightest chargino as 
a function of the mass for tan /3 = 10 and /i > 0. 

Clearly, the invisible decay of the neutralino into neutrinos cannot be reconstructed and for 
this reason we discarded it. On the other hand, it is possible to measure the chargino decay 
vertex in the decays: 

• Xi ~^ ^q'q. denoted jj; 

• xt ^ T^QQ^ called rjj; 

• xf ^ T^bb, called rbb; 

• Xt ^ r+£+£-, called tU; 

• xt ^ (^~^bb, that we denote by £bb; 

• xt ^ that we denote by iii; 

• xt ^ that we denote by ijj. 

Although the decay Xi ~* vv£^ can give rise to a or //^ with a high impact parameter, 
we do not consider the mode in our study since it is not possible to obtain its decay vertex. 

We considered a crude model of the LHC detectors in order to identify events containing 
displaced vertices. We selected events presenting neutralino or chargino decays away from 
the primary vertex via the requirement that the displaced neutralino/chargino vertex is 
outside an eUipsoid around the primary vertex 
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where the 2;-axis is along the beam direction. We took 5xy = 100 /xm and 5z = 2.5 mm 
that correspond to five times the expected resolution in each direction. To guarantee a high 
efficiency in the reconstruction of the displaced vertices without a full detector simulation, 
we required the tracks leaving the displaced vertex to be inside the rapidity coverage of the 
vertex detector, i.e. \ri\ < 2.5. Moreover, we also required that the displaced vertices are 
inside the vertex detector - that is, the vertices must be within a radius of 550 mm and 
z-axis length of 800 mm. 

The SM backgrounds coming, for instance, from displaced vertices associated to fe's or 
r's can be eliminated by requiring that the tracks defining a displaced vertex should have 
an invariant mass larger than 20 GeV. This way the displaced vertex signal passing all the 
above cuts is essentially physics background free, however, there might exist instrumental 
backgrounds which were not considered here. 

An important issue in the displaced vertex search is the trigger on the events containing 
them. In order to mimic the triggers used by the LHC collaborations, we accept events 
passing at least one of the following requirements: 

• The event has one isolated electron with px > 20 GeV; 

• The event has one isolated muon with > Q GeV; 

• The event has two isolated electrons with pr > 15 GeV; 

• The event has one jet with pt > 100 GeV; 

• The event has missing transversal energy > 100 GeV. 

IV. DISCOVERY REACH AT THE LHC 

We estimated the LHC discovery reach in all the channels described in Sec. IIII A[ We 
required that either the signal leads to 5a departure from the background where this is not 
vanishing or 5 events in regions where there is no SM background. We present our results 
in the ® plane for tan/3 = 10, > and integrated luminosities of 10 fb~^ and 100 
fb-i. 

We depict in Figure [9] the LHC discovery potential in the all inclusive channel (IN) . For 
the sake of comparison, we also present in this figure the discovery reach assuming R-parity 
conservation. As we can see from this figure, the introduction of bilinear R-parity violation 
reduces the reach in 7723/2 for a given value of mo if we use the search strategy design for the 
R-parity conserving case. Basically, the decay of the LSP reduces the missing transverse 
energy making it harder to disentangle the SUSY signal and the SM background. Comparing 
the left and right panels of Fig. [9], we can see that a larger luminosity in this channel expands 
the reach in mg/a from ~ 70-80 TeV to ~ 90-100 TeV. 

The reach in the li channel is presented in Figure [TOl In the R-parity conserving scenario 
this channel is sensitive to the spectrum details of mAMSB since the lightest chargino 
decay contains soft charged particles in addition to the LSP and consequently, there are 



19 




500 1000 1500 2000 2500 3000 3500 4000 500 1000 1500 2000 2500 3000 3500 



m„(GeV) m„{GeV) 

FIG. 9: Discovery reach in the inclusive channel (IN) for tan/3 = 10, fi > and the integrated 
luminosities of 10 fb~^ (left panel) and 100 fb~^ (right panel) . The dark (blue) square mark the 
points where niAMSB with R-parity conservation can be discovered while the stars stand for the 
reach of our BRpV-mAMSB model. Dark (red) circles stand for the points excluded by LEP data 
while no solution for the neutrino masses was found in the grey (cyan) squares. The dark (red) 
area on the left is theoretically excluded by the existence of tachyonic particles and the light shaded 
area adjacent to this one is the area with a stau LSP. The light (green) line is the contour for the 
gluino mass of 1 TeV. The Higgs mass is 114 GeV on the dashed dark line. 



less channels available that give rise to hard leptons. Once again the introduction of R- 
parity violation depletes this signal because of the reduced missing transverse momentum. 
Moreover, the extra produced leptons in the chargino or neutralino decays have the tendency 
of contributing to the trilepton topology. Comparing the left ad right panels of this figure 
it is clear that the reach in this channel is extended with the increase of the luminosity, 
however, its reach is still smaller than the IN one. 

The reach in the multilepton channel is depicted in Figure [TTl For 10 fb~^ the R-parity 
conserving scenario has a very limited reach with the signal being sizeable only in the area 
presenting light sleptons. The inclusion of BRpV interactions increases the LHC reach in 
the M£ channel at small values of with the IN and M£ having similar reaches at low 
luminosities. At higher luminosities, 100 fb~^, the multilepton channel reach is considerably 
extended with and without R-parity conservation, being this the channel with largest reach 
at small mo. Moreover, this is the SUSY canonical topology with the largest potential for 
discovery in the R-parity violating scenario. 

We also studied the exclusive channels containing two isolated charged leptons verifying 
that the introduction of R-parity violating interactions enhances these channels. In the 
R-parity violating scenario the LHC OS reach is similar to the 1£ channel one except at 
large mo where the OS signal has a reduced reach; see Figure [121 On the other hand, the 
SS channel presents a smaller SM background in addition to an enhanced signal due to the 
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Discovery reach in the multilepton channel for the parameters and conventions used in 



presence of majorana states in SUSY models. Therefore, this topology has a good discovery 
reach in our BRpV-mAMSB model as can be seen from Figure [T3l In fact, the SS final state 
has a slightly larger reach than the fully inclusive mode IN, with the SS channel being the 
second most important channel in the presence of BRpV interactions. 

Before we move on to the displaced vertex signal, we would like to stress that our re- 
sults are an indication of R-parity violating interactions in the canonical SUSY searches. 
Certainly, the introduction of a larger number of floating cuts leads to larger reaches in all 
channels, like the ones in Ref. 29|]. Nevertheless, the depletion of the fully inclusive channel 
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FIG. 12: Discovery reach in the opposite sign lepton channel for the parameters and conventions 
used in Fig. [H 

and the enhancement of exclusive topologies containing isolated leptons, as compared to the 
R-conserving case, must persist in a more elaborate analysis. In this light, it should be 
possible to determine whether R parity is broken or not by combining the results of all the 
canonical SUSY search channels. 
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FIG. 13: Discovery reach in the same sign lepton channel for the parameters and conventions used 
in Fig. El 



The smoking gun of BRpV-mAMSB is the existence of detached vertices exhibiting a high 
invariant mass. There is no SM background for these events except for possible instrumental 
backgrounds, rendering this channel a very strong evidence for Physics beyond the standard 
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model. We present in Figure [T3] the reach in the displaced vertex topology for integrated 
luminosities of 10 and 100 fb~^. As we can see from this figure, this channel does have the 
largest reach 110 (120) TeV for 10 (100) fb~^) with the nice feature of being almost 
independent of mo, except at small mo where the presence of light scalars lead to rapid 
neutralino and chargino decay; see Fig. [HI 
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FIG. 14: Discovery reach for the LHC in displaced vertices channel for tan /3 = 10 and ^Li > in 
the ® mo plane and an integrated luminosities of 10 fb^^ (stars) and 100 fb^^ (dark (blue) 
squares). The light (green) lines gives the contour of a gluino mass of 1 TeV and 2 TeV. All other 
conventions follow the ones of FiglOl 



V. CONCLUSIONS 

We have studied the phenomenology of AMSB model augmented with bilinear R-parity 
parameters at the LHC. We show that the presence of bilinear R-parity interactions modifies 
the canonical channels used for the supersymmetry search. The decay of the neutralino and 
chargino weaken the fully inclusive signal, however, the existence of further leptons in the 
final state leads to an enhancement of the leptonic exclusive modes. In the case that the 
BRpV-mAMSB final state contains three or more charged leptons the reach is so enhanced 
that this channel alone has a comparable reach to the fully inclusive one with R-parity 
conservation. One interesting aspect of the drastic change in the reach of the different 
topologies is that a possible positive signal at the LHC can be used to disentangle models 
with and without R-parity conservation. 
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FIG. 15: Discovery reach at the Tevatron in displaced vertices channel for tan /3 = 10,/i>0in the 
"^3/2 ^ "^0 plane and an integrated luminosity of 2 fb^^ (stars) and 10 fb~^ (dark (blue) squares). 
The light (green) lines gives the contour of a gluino mass of 1 TeV and 2 TeV. All other conventions 
follow the ones of FigjUJ 

Our BRpV-mAMSB model connects the R-parity violating parameters to measured neu- 
trino properties. Due to the smallness of the couplings needed to reproduce the observed 
neutrino masses and mixings, our model predicts a rather large lifetime for charginos and 
neutralinos. This leads to the smoking gun of the BRpV-mAMSB theory, that is, the exis- 
tence of displaced vertices associated to neutralino or chargino decays that exhibit a large 
visible invariant mass. We showed that the search for such detached vertices does lead to 
the largest discovery reach for such models. Indeed this channel can probe m3/2 up to ^ 100 
TeV, which constitutes almost all the natural parameter space. At this point it is interesting 
to point out that the search of detached vertices at the Tevatron can cover a large fraction of 
the presently allowed parameter space. We depict in Figure [15] the parameter space region 
that can be probed at the Tevatron for an integrated luminosity of 2 fb~^, showing that the 
Tevatron collaborations can already cover a lot of ground in searching for BRpV-mAMSB. 
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